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We propose a dynamical explanation of the Meissner effect in superconductors and predict the 
existence of a spin Meissner effect: that a macroscopic spin current flows within a London penetration 
depth Al of the surface of superconductors in the absence of applied external fields, with carrier 
density = the superfluid density and carrier speed v = ft/(4m e Ai) (m e =bare electron mass). 
The two members of a Cooper pair circulate in orbits of radius 2Al in opposite direction and the 
spin current in a Cooper pair has orbital angular momentum h. Our description also provides a 
'geometric' interpretation of the difference between type I and type II superconductors. 



PACS numbers: 

What is the force that causes the electrons near the 
surface of a metal to start moving and develop the electric 
current that will screen the magnetic field in the interior 
when the metal is cooled into the superconducting state 
in the presence of an external static magnetic field? De- 
spite having been known for over 70 years, no dynamical 
explanation of the Meissner effect exists within the con- 
ventional understanding of superconductivity [l|. 

Superfluid electrons carry one Bohr magneton of mag- 
netic moment in addition to one electron charge. One 
may similarly expect that the superfluid will respond to 
an electric field by developing a spin current^,^. We 
assume there is no externally applied electric field, never- 
theless superfluid electrons interact with the electric field 
originating in the positive ionic medium where they re- 
side. Spin currents do not break time reversal invariance 
and hence should be insensitive to non-magnetic disorder. 
Will superconducting electrons respond to the presence 
of internal electric fields by developing a spin current? 

In this paper we propose answers to both of these 
questions and in addition we offer a new interpreta- 
tion of what a Cooper pair is. We rely on key physics 
proposed within the theory of hole superconductivity 
that superconductivity is driven by kinetic energy low- 
ering, that is achieved by expansion of the electronic 
wave function [f|, that also results in expulsion of nega- 
tive charge from the interior towards the surface[6[. 

Consider a metal in a uniform magnetic field B along 
the z direction. Its diamagnetic response (Landau dia- 
magnetism) can be understood as arising from induced 
Ampere circular currents of radius r ~ kp l . For a nearly 
filled band (as required in the theory of hole supercon- 
ductivity), r ~ a, the lattice spacing[5]. Assume the 
transition to superconductivity involves a radially out- 
ward motion of the electrons in a Cooper pair to a radius 
r = 2Xl. The electrons will acquire an azimuthal veloc- 
ity v<p as they move outward due to the magnetic Lorentz 
force Fl = (e/c)v x -B[l[. Conservation of angular mo- 
mentum dictates that the azimuthal velocity at radius 
r is given by: 
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and for radius r = 2Al 

V4, 



e\L - 
1 

m e c 



(2) 



Note that the magnetic vector potential at position r 
for a uniform magnetic field is given by A = B x r/2, 
hence Eq. (1) can be written as 



V& = A 

m e c 



(3) 



In the conventional theory, Eq. (3) is understood as aris- 
ing from the 'rigidity' of the superconducting wavefunc- 
tion, which causes the canonical momentum p in the ex- 
pression 



m e v = p A 

c 



(4) 



to remain independent of the value of the applied mag- 
netic field Q. Here instead, we derived this result from 
the dynamics. 

When we now consider all the Cooper pairs in the su- 
perconductor, the 'internal' azimuthal velocities get can- 
celled by superposition, as shown in Fig. 1(a), however 
those near the surface are not. The uncompensated az- 
imuthal velocity Eq. (2) of the carriers near the surface 
is precisely what is needed to screen the magnetic field 
in the interior of the superconductor. For a long cylinder 
with superfluid density n s the current density is 



J charge — J en s V§. 



(5) 



Such current density flowing in a surface layer of thick- 
ness Al induces a magnetic field in the interior of mag- 
nitude 
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and direction opposite to B, that will exactly cancel B 
since H 
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FIG. 1: Illustration of the Meissner and spin Meissner effects, 
(a) shows the net direction of orbital motion of the negatively 
charged members of the Cooper pair that expanded to radius 
2Al (small circles in (a)), arising from the Lorentz force due 
to a magnetic field pointing out of the paper. Note how the 
'internal' velocities get cancelled and there is a net counter- 
clockwise motion near the surface, (b) shows the spin Meiss- 
ner effect for one of the circles in (a) : up and down magnetic 
moments get deflected in opposite direction due to the torque 
exerted by the radially pointing electric field from the posi- 
tive charge distribution, (b) also shows the equivalent electric 
dipoles as they move out. Internal dipoles cancel out and only 
the dipoles near the surface of the sample survive. 



Having thus provided a 'dynamical' explanation of the 
origin of the ordinary Meissner effect, we now turn our at- 
tention to the spin Meissner effect. We choose the z axis 
as the spin quantization axis, hence the electron magnetic 
moment is 



p = ps&z 
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with a = +/ — I. Consider now the effect of a uniform 
positive background charge density on a radially outgo- 
ing magnetic moment p (Fig. I (b)). The moving mag- 
netic moment is equivalent to an electric dipole p in the 
laboratory frame [£j, [lfj 



p = 7— x p 
c 



(9) 



(7 ~ 1) on which a radial electric field Ei will exert a 
torque f = p x Ei, causing a change in orbital angular 
momentum 

dL d ,v . t> , . 

—r = m e — (r x v ) — - x p x E t 10 
dt dV ' c p; v ' 

The electric field of a cylindrical uniform charge density 
p is Ei — 2npr. The positive ionic charge density that 
the superfluid senses is p = \e\n s . Hence from Eq. (10) 

_ dv ,v . . 

m e r x — = 2-Ken 8 r X (— X p) (11) 
dt c 

Eq. (11) can be interpreted as if an effective spatially 
uniform mag netic field|l3, EH 



B„ 



2-Kn s p 



(12) 



is exerting a Lorentz force F = (e/c)v x on the elec- 
tron. Just as in the earlier discussion, it will give rise to 
an azimuthal velocity 



V 0a 



-B„ 
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for outward motion to radius r — 2Al, where now the 
direction of the azimuthal velocity is opposite for up and 
down spin electrons. Again, the 'internal' azimuthal ve- 
locities get cancelled by superposition and the uncom- 
pensated azimuthal velocities give rise to a macroscopic 
spin current 



J spin 



y(w 0T - Vol) = n sV01 
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flowing within a London penetration depth of the surface, 
whether or not an external magnetic field is applied. 

Just like in the case of the real magnetic field, the 
effective magnetic field Eq. (12) can be understood as 
arising from a vector potential 



A„ 



B a xr 
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and the total velocity 



v a = A A a = vs + v 0o 

m e c m e c 



(15) 
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can be interpreted as arising from the 'rigidity' of the 
wave function as in Eq. (4). 

Note that the Dirac Hamiltonian for an electron in- 
cluding the spin-orbit interaction in the non-relativistic 
limit in the presence of an electric field E yields [l2| 
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The last term in Eq. (17) has been interpreted as a mag- 
netic vector potential giving rise to an 'effective' magnetic 
field. [Til. [131 However, with E = Ei it yields an effective 
magnetic field of opposite sign to Eq. (12), so this argu- 
ment is misleading. That the force on an electron with 
magnetic moment p moving with velocity v in a positive 
background is in direction ev x p, as predicted by Eq. 
(12), also follows from the treatment of Ref. [Ioj |. 

Using Eq. (7), the spin current velocity Eq. (13) can 
be written in the remarkably simple form 



V a 



4m P \] 



(18) 



which is the central result of this paper. 

How large is this azimuthal velocity? For a typical 
X L = 4001, Eq. (18) yields v 0a = 72,395cm/ s. The 
kinetic energy per carrier associated with it is 
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which yields lA9/ieV for Al = 40(M, which is compa- 
rable to superconducting condensation energies per elec- 
tron. Since each superfluid electron in the bulk lowers 
its energy by the condensation energy, there is plenty 
of energy available to account for the kinetic energy of 
the spin current that resides only in the surface layer of 
thickness A^. Note also the form of the last term in Eq. 
(18) (obtained using Eq. (7)): it says that for two oppo- 
sitely alligned spin 1/2 magnetic moments perpendicular 
to the radius joining them at distance of order the inter- 
electron distance (nj 3 ), the energy lowering from mag- 
netic dipole-dipole interaction can account for the spin 
current kinetic energy cost. 

There is another independent argument that supports 
the validity of Eq. (18). An applied magnetic field B = 
Bz will exert a Lorentz force on the superfluid carriers 
in the spin current in the surface layer. This force 



Ft. 



-v Qa x B 
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(20) 



points inward for magnetic moment parallel to B (a = 
— 1) and outward for opposite spin. In addition, because 
B only penetrates a distance ~ 2Al[I3|, its gradient ex- 
erts a force on the magnetic moments that is outward for 
magnetic moment parallel to B and inward for opposite 
spin: 



F„ 



eh dB „ 
2m e c dr 
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where we have used dB/dr ~ Bj 2A^. If we require that 
the forces Eq. (20) and (21) exactly cancel for each spin 
component we obtain the value Eq. (18) for the spin 
current velocity. Hence an applied B will not affect the 
radial distribution of the spin current. Note how in this 
argument the orbital and spin aspects are inextricably 
intertwined. 

What is the magnitude of a real magnetic field that 
would give rise to the same azimuthal speed of the su- 
perfluid as the speed of the spin current? Equating Eqs. 
(2) and (18) it follows it is simply \B a \ (Eq. (12)), which 
can also be written as 



\BJ = - 



he 
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with 4>q — hc/2\e\ the flux quantum. Application of an 
external magnetic field of magnitude B s would cause one 
of the components of the spin current to move twice as 
fast and the other one to stop altogether. However, note 
that this magnetic field is essentially the lower critical 
field of a type II superconductor [§[ 



Inn 



(23) 



with k the Ginzburg Landau parameter (k > l/v2 for 
type II superconductors). Hence, when a magnetic field 
is applied to a type II superconductor large enough that it 





(a) Type I superconductor 



(b) Type II superconductor 



FIG. 2: In type II materials ((b)) with £ < 2Al, a normal 
vortex core of diameter £ can be enclosed by both orbits of 
the same Cooper pair. In type I materials ((a)), £ > 2A_l and 
this cannot happen. 



causes one of the spin current components to essentially 
stop (B — B s ~ H c i), superconductivity is destroyed 
and the magnetic field penetrates the sample forming flux 
tubes of normal material threaded by magnetic field H c \ . 

What about type I superconductors? The thermody- 
namic critical field for a free electron model can be writ- 
ten as 8] 



Br 



3 1 



2 7T 2 £oA L 



(24) 



with £o the Pippard coherence length. Hence from Eqs. 
(22) and (24), B s = 0.64(£ /A L )# C . Since for a type I 
superconductor £o > V2Xl, the superconducting state 
will be destroyed well before the applied magnetic field 
reaches B s . 

What is the role of the coherence length £? We pro- 
pose that £ is the average distance between the centers of 
the orbits of the two members of a Cooper pair, as shown 
schematically in Fig. 2. This is consistent with the con- 
ventional understanding that £ is the 'size' of the Cooper 
pair. However, in our case it provides a new geomet- 
ric interpretation for why the cross-over between type I 
and type II behavior occurs for £ ~ 2A^ . This will be 
discussed in more detail in a forthcoming paper. 

How can one detect the existence of this spin current 
in superconductors experimentally? (i) Inserting in a su- 
perconducting ring a spintronics device that rectifies spin 
current [15j one would measure an unexpected asymmetry 
depending on its orientation; (ii) Spin-resolved neutron 
scattering with very cold neutrons may be able to de- 
tect it; (iii) Since the spin current gives rise to equivalent 
electric dipoles, a superconducting ring will have electric 
field lines as shown in Fig. 3, with a non-zero quadrupole 
moment which we can estimate as 



Q = —4:7ra 2 h\i J n s p 



(25) 



with a the radius, h the height, n s the superfluid density 
and p the electric dipole moment Eq. (9). For Ax = 
400A, p = 2.24 x 10~ 26 esu, n s = 1.77 x 10 22 /cm 3 and 



Q 



-1.99 x 10 s a(cm) 2 h(cm)esu 



(26) 



One may attempt to measure the electric field around 
the ring or the electrostatic force between two such rings. 
For two rings on top of each other (Fig. 4a) a distance 
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FIG. 3: Pattern of electric field (curved lines) expected to 
arise for a superconducting ring. The vertical arrows denote 
the direction of the magnetic moment of the superelectrons 
and the horizontal arrows the direction of their motion. No 
magnetic field nor charge current is present. 




FIG. 4: Force between two superconducting rings due to the 
electric fields produced by the spin currents. Depending on 
their relative position and orientation the force can be repul- 
sive (a) or attractive (b), (c). In (d), a smaller ring should be 
in stable equilibrium at the center of a larger ring. 

z apart the force is repulsive and given by F z — 5Q 2 /z 6 
and should be detectable for small enough z (note that 
z > h). For different relative position and orientation 
of the rings the force will be repulsive or attractive as 
shown in Fig. 4. 

In a spherical or planar geometry in the absence of 
external fields there is no preferred direction to the spin 



current. However applying a small magnetic field will 
define a spin quantization direction and orient the spin 
current flow according to Eq. (14). An externally applied 
electric field is expected to modify the spin current at low 
temperatures. We point out that the magnitude of the 
spin current density ej sp i n predicted by this physics can 
be as high as 10 8 Amps /cm 2 and may provide a useful 
source of spin current for spintronics applications. 

Why do supcrfhiid carriers near the surface 'choose' to 
pay the kinetic energy Eq. (19) in the absence of exter- 
nally applied fields, rather than stop and save that en- 
ergy altogether? Those carriers are tasked with shielding 
the interior of the superconductor from applied magnetic 
fields by setting up a Meissner current. Thus we may 
think of them as 'confined' in the surface layer of thick- 
ness ~ 2\l, ready to respond. The quantum zero-point 
motion energy of an electron confined in a region of lin- 
ear dimension Ax is h 2 /2m e (Ax) 2 , of the same order as 
the kinetic energy Eq. (19). 

Why is the radius of the orbit 2Al? 2Al = \Jhcj\cB a \, 
the spin-orbit 'magnetic length'. The flux of B a through 
the circle of radius 2A^ is 




the flux quantum. The spin current velocity Eq. (18) cor- 
responds to an orbital angular momentum m e v 0a r = ft/2 
for an electron in a circular orbit of radius r = 2A^. As 
the Cooper pair is born, its wave function expands so 
that the orbit of each member of the Cooper pair will 
thread one flux quantum of spin-orbit interaction flux, 
and in the process each member of the Cooper pair ac- 
quires orbital angular momentum h/2, equal to its spin 
angular momentum. Thus the Cooper pair wave func- 
tion is single-valued, and the factor of 2 in (j) = hc/2e 
acquires a new meaning. 
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